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Abstract. The phenomenon of evanescent wave amplification inside of a half- 
wavelength-thick wire medium slab is studied numerically and used for subwavelength 
imaging of objects at significant distances from the lens in this paper. The wire 
medium is analyzed using both a spatially dispersive finite- difference time-domain 
(FDTD) method employing the effective medium theory and a full-wave commercial 
electromagnetic simulator CST Microwave Studio™ modelling the physical structure. 
It is demonstrated that subwavelength details of a source placed at a distance of 
A/ 10 from a wire medium slab can be detected inside the slab with a resolution of 
approximately A/ 10 in spite of the fact that they cannot be resolved at the front 
interface of the device, due to the rapid decay of evanescent spatial harmonics in free 
space. The influence of different distances between the source and the wire medium 
on the internal imaging property is also addressed. 
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1. Introduction 

Conventional imaging systems are restricted by the so-called diffraction limit: any source 
details below the half-wavelength at the frequency of operation cannot be resolved at 
the image plane. However, it has been proposed in [1] that a planar lens formed by 
the left-handed material (LHM) [2] can be used to image source information with 
a spatial resolution below the diffraction limit. A lens formed by the LHM with 
s = /i = — 1 is termed as the 'perfect lens' and the principle of its operation is based on 
the negative refraction of propagating waves and the amplification of evanescent field 
components [1]. It is well known that evanescent waves carry subwavelength source 
details and decay exponentially in positive-index materials. The uniqueness of the 
LHM is that the evanescent waves show a growing behavior when propagating within 
the medium. Hence the subwavelength details which are lost in the region between the 
source and a perfect lens can be restored to create a perfect image [1]. The amplification 
of evanescent waves in the LHM is due to the resonant excitation of surface plasmons 
at the interfaces. However, such effect is sensitive to the losses in the LHM thus limits 
the maximum thickness of the LHM slab [3]. Furthermore, the mismatch of the LHM 
with its surrounding medium also limits the imaging capability of LHM lenses [4] . 

Recently it has been suggested to use an alternative way to transport subwavelength 
source details to an image plane at a significant distance. The principle of operation, 
referred as 'canalisation' [5], is based on that for a certain type of devices, the evanescent 
wave components can be transformed into propagating waves, and therefore the source 
field can be delivered to its back interface with little or no deterioration. In contrast 
to the case of the LHM, such devices are less sensitive to losses. Recent works on 
the subwavelength imaging using anisotropic materials (operating in the canalisation 
regime) include [6-8]. Moreover, the work presented in [9] shows that introducing 
nonlinearity is able to suppress the diffraction, and is useful for future subwavelength 
imaging devices even when the dispersion curve is not flat. 

One typical example of the structures operating in the canalisation regime is the 
wire medium formed by an array of parallel conducting wires [10]. The thickness of the 
wire medium needs to be equal to an integer multiple of half- wavelengths at the operating 
frequency (due to the Fabry-Perot resonance) in order to avoid reflections between the 
source and the structure. It has been demonstrated experimentally and numerically 
that such a canalisation regime does indeed exist and, subwavelength details can be 
transported through the wire medium. 

In addition to the canalisation regime, another interesting phenomenon of the 
wire medium is observed from numerical simulations and presented in this paper: the 
evanescent waves can be amplified inside the wire medium. Hence the effect of images 
appearing inside of the wire medium slab is called 'internal imaging'. The previously 
investigated subwavelength imaging property of the wire medium, i.e. images appear at 
the back interface, is highly dependent on the distance between the source and the wire 
medium, and the external images disappear when the distance increases. Due to the 
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internal imaging effect, the subwavelength imaging capability of the wire medium can 
be further improved. In contrast to the case of LHM lenses where the amplification of 
evanescent waves is due to the resonant excitation of surface plasmons, the amplification 
of evanescent waves in a wire medium slab is a result of the resonant excitation of 
standing waves inside the slab. Although the internal imaging effect also exists in LHM 
slabs (there are two focal points at both inside and outside of LHM lenses), as mentioned 
above, the practical realization of LHM lenses has many difficulties. On the other hand, 
the internal imaging capability of the wire medium can be used for subwavelength 
imaging of objects at significant distances from the lens. In the following, we first 
present the numerical simulation results using a spatially dispersive finite- difference 
time-domain (FDTD) method developed in [11], where the effective medium theory 
is used and micro-structures of the wire medium are neglected. Hence the FDTD 
method has better efficiency comparing with simulations modelling micro-structures. 
Then the FDTD simulation results are validated by a commercial simulation package 
CST Microwave Studio™ modelling the actual physical structure of the wire medium. 



2. Spatially dispersive FDTD simulations of internal imaging by the wire 
medium 

The simulation of the wire medium can be performed either by modelling the physical 
structure i.e. parallel conducting wires, or using the effective medium theory, so that 
the wire medium can be represented by a homogeneous dielectric slab provided that the 
inner spacing of wires is much smaller than the wavelength of interest. In this section, we 
demonstrate the internal imaging phenomenon numerically using the effective medium 
theory by applying the spatially dispersive FDTD method [11]. modelling the physical 
structure based on a commercial simulation package CST Microwave Studio™ is 
introduced in the next section. 

The wire medium can be approximated as a homogeneous dielectric slab with both 
frequency and spatial dispersions. In this paper we use the following expression for the 
permittivity tensor derived in [12]: 

k 2 

e(u, q) = xx + e(u, q y )yy + zz, e(u, q y ) = 1 - p (1) 

K -q y 

where the y-axis is oriented along wires, q y is the y-component of the wave vector q, 
k = uj/c is the wave number of the free space, c is the speed of light, and k p is the wave 
number corresponding to the plasma frequency of the wire medium. The wave number 
k p depends on the lattice periods a and 6, and on the radius of wires r [13]: 

p l„g + F( /6)' K ' 2 tA n ) 6 ? < ' 
For the case of the square grid (a = 6), F{X) = 0.5275 and the expression (ED reduces to 
2 = 2n/a 2 

v i n a +Q.5275' V ' 

ZTTT 
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The expression (OQ) for the permittivity tensor of the wire medium is valid, and the 
structure is non-magnetic if the wires are thin as compared to the lattice periods (when 
the polarization across the wires is negligibly small as compared to the longitudinal 
polarization) and if the lattice periods are much smaller than the wavelength (when the 
wire medium can be homogenized). 

A spatially dispersive FDTD method has been developed [11] to take into account 
the spatial dispersion effect in the wire medium [12]. We have applied the spatially 
dispersive FDTD method to analyze the subwavelength imaging property of the wire 
medium [11]. In this paper, we extend the method to three dimensions to study the 
internal imaging property of the wire medium. 

The computation domain in FDTD is shown in Fig. [1] (except that the wire medium 
is modelled as a bulk material instead of the parallel thin wires). A crown-shaped near- 
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Figure 1. A three-dimensional (3-D) FDTD simulation domain for the modelling of 
internal subwavelength imaging by a wire medium slab with thickness A/2 and infinite 
transverse dimensions. 

field source is applied in the simulations and it is placed at a distance of d to the 
front interface of the wire medium. The distance d varies from A/15, A/10 to A/7. 
The source is modelled as segments of perfect electric conductor (PEC) using staircase 
approximations, and it is excited at the bottom left corner using a delta-gap excitation 
in FDTD [14]. The radius of the PEC wires is the same as the FDTD cell size, which 
is fixed at A/ 100 for all FDTD simulations. The time step is chosen according to 
the Courant stability criterion [14]. Berenger's perfectly matched layer (PML) [15] is 
employed to truncate the outer boundary of the simulation domain along ^-direction 
and the modified PML [16] is applied to model infinite wire medium slabs along x- and 
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^-directions and reduce the convergence time in simulations [16]. The dimensions of 
the simulation domain in x- and z-directions are 2 A/3 x 2 A/3. The thickness of the 
wire medium slab is chosen as A/2 to achieve good impedance matching between the 
source and the structure. As shown in Fig. [U the middle plane of the wire medium is 
defined as the internal image plane. It has been demonstrated that the wire medium 
can image the electric field component along the wires [10]. Inside the wire medium, the 
field distribution has TEM polarization (transverse electric and magnetic modes with 
respect to the orientation of wires) and it is formed by the so-called transmission-line 
modes [5]. The electric field components of such modes are zero along the wires (but 
the electric displacement is non-zero). Therefore, in Fig. [2] we plot the distributions of 
electric field only at the source plane and the front interface, but at the internal image 
plane the electric displacement (normalized to the free space permittivity) is plotted 
instead. It is equivalent to the electric current density flowing along the wire since 
D y = £oE y + J y /(juj). The field dist r ibut ions at various planes are calculated according 
to different distances d between the source and the wire medium and plotted after the 
steady state is reached in simulations. 
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Figure 2. Normalized distributions of the ^-component of the electric field at the 
source plane, the front interface, and the electric displacement at the internal image 
plane of the wire medium for different distances between the source and the wire 
medium calculated from spatially dispersive FDTD simulations. 



It is clearly seen that the crown-shaped source creates a subwavelength field 
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distribution at the source plane, however, it becomes blurred at the front interface 
of the wire medium due to the spatial decay of the evanescent spatial harmonics in the 
free space region, especially for the case of d = A/7 that the distribution at the front 
interface can be hardly identified. The wire medium operates in the canalisation regime 
and the electric field distribution at the back interface is repeated at the front one, 
however, such field distribution does not create a sharp image of the original source. 
The surprising result is that inside the wire medium, lost subwavelength details can be 
recovered due to the amplification of evanescent wave components. In contrast to the 
case of LHM slabs, all incident evanescent spatial harmonics inside the wire medium are 
transformed into propagating waves, independent on their transverse variations. The 
thickness of the wire medium slab is tuned to the Fabry- Perot resonance. As a result, 
the evanescent waves excite standing waves inside the slab and in this case, they are 
amplified due to the natural resonance of the wire medium slab. It can be identified 
from Fig. [2] that the resolution of the internal images is about A/10. 

The internal images with good subwavelength resolutions are observed when the 
source is placed close to the wire medium slab (i.e. d < A/ 15). If the distance is increased 
to A/10 and even A/7, the internal images do not contain enough subwavelength details 
of the source. We conclude that internal imaging with good subwavelength resolution is 
available only for a limited range of distances between the source and the wire medium, 
because the dependence of the transfer function for the evanescent spatial harmonics on 
the transverse wave vector is linear rather than exponential as for the case of LHM [1]. 

In the above analysis, we have assumed the wire medium to be a homogeneous 
dielectric material. However in practice, since the electric displacement inside the wire 
medium is proportional to the currents flowing along the wires, in order to detect the 
internal images, one needs to measure the current density instead of either electric or 
magnetic field, in contrast to the canalisation regime that the image is formed at the back 
interface of the wire medium and can be detected using the near-field scanning. In the 
following section, we model the physical structure of the wire medium and demonstrate 
numerically that the currents along the wires indeed carry the subwavelength details of 
a source. 

3. Simulations of internal imaging by modelling the physical structure of 
the wire medium 

In the simulations using CST Microwave Studio™, the overall dimensions of the wire 
medium and the crown-shaped source are kept the same as in the previous section 
using the FDTD method. The number of wires is 21 x 21 along x- and z-directions 
respectively, and the inner spacing between wires is a = A/ 30. The radius of the wires 
is r = A/ 300. We also consider three different distances between the source and the 
wire medium: d — A/15, d — A/10 and d = A/7. The normalized distributions of the 
^/-component of the electric field at the source plane, at the front interface, as well as 
the current density along the wires (the currents are converted to field intensities to 
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compare with the other two distributions) in the middle plane of the wire medium are 
shown in Fig. [31 It is shown that the source distribution is much sharper compared to 
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Figure 3. Normalized distributions of the ^-component of the electric field at the 
source plane, the front interface, and the current density at the internal image plane 
of the wire medium for different distances between the source and the wire medium 
calculated from CST simulations of the physical structure. 



the FDTD simulation results. This is due to the fine mesh used in CST simulations 
in order to model thin wires accurately. Again we observe similar distributions at the 
front interface of the wire medium where less details of the source can be detected. The 
current density is plotted on the discrete grids (the third column in Fig. [3]) since only 
21 x 21 wires are present in the modelled wire medium. The current density shows 
significant enhancement compared to the field distributions at the front interface of the 
wire medium, which demonstrates the internal subwavelength imaging capability of the 
wire medium. 

The low image resolution of the internal current density distribution is due to 
the insufficient sampling points (21 x 21). Therefore to acquire a higher resolution of 
internal images, one can use the wire medium formed by more densely packed wires. For 
example, we have considered a wire medium formed by 41 x 41 wires with the spacing 
between the wires a = A/60, keeping the transverse dimension of the wire medium the 
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same as the previous case. According to Eq. ([2j), the plasma frequency is higher for 
this case and therefore one could expect images with a higher resolution. However the 
wire density cannot be increased arbitrarily due to the fact that when the filling ratio of 
wires increases, the structure may exhibit anisotropic magnetic properties [17-19]. The 
influence of the material parameter change on the subwavelength imaging property of 
the wire medium is currently under our investigation. 

The current density distributions from CST simulations for the case of 41 x 41 
wires are plotted in Fig. 0J Note that for illustration purposes, the distributions are not 
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Figure 4. The current density from CST simulations of the physical wire medium 
(41 x 41 wires) for different distances between the source and the wire medium. 

normalized. It can be seen that by increasing the density of wires, the internal image 
resolution can be improved, which is especially useful when a more complex source is 
considered. 

The preliminary estimations show that the image appears near the middle plane of 
the A/2 thick wire medium slab and does not change significantly even when the image 
plane is offset from the middle plane by 5-10%. The plane of the best image formation 
cannot be determined precisely from numerical simulations since the distribution of 
electric displacement or current density inside the wire medium varies slowly along the 
direction of wires. We have considered a wire medium slab with a half-wavelength 
thickness in this paper. For the wire medium slabs with greater thicknesses (an integer 
multiple of a half- wavelength) , there are multiple planes for the internal imaging, where 
the evanescent spatial harmonics generated by the source can be amplified and the 
information lost in the free space region can be restored. The number of internal image 
planes can be given by TV if the thickness of the wire medium is NX/2. The distance of 
the internal image planes to the front interface of the wire medium can be approximately 
given by (2n — l)A/4 where n — 1, 2, N. 

4. Conclusion 

The internal imaging capability of the wire medium, i.e. images with a much 
higher resolution appear inside of the wire medium comparing with the external 
ones, is demonstrated numerically by modelling the physical structure and using the 
effective medium theory: the wire medium slabs are capable of imaging a source with 
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subwavelength resolutions located at significant distances away from the device. The 
operational principle of 'internal imaging' is based on the excitation of standing waves 
inside the wire medium slab rather than the resonance of surface plasmons in LHMs. The 
features of being low loss and easy in manufacturing make the wire medium attractive 
for practical applications. In practice, it is possible to capture the internal image 
by embedding detectors into metallic wires at the internal image plane and directly 
measuring the currents flowing along the wires. The internal imaging effect can be 
used for subsurface imaging and the creation of new generation subwavelength imaging 
devices. 
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